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Abstract 

Unavoidable variations in size and position of the building blocks of photonic crystals cause 
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Q,! light scattering and extinction of coherent beams. We present a new model for both 2 and 3- 
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• ' dimensional photonic crystals that relates the extinction length to the magnitude of the variations. 
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•'-H ■ The predicted lengths agree well with our new experiments on high-quality opals and inverse opals, 
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(— I , and with literature data analyzed by us. As a result, control over photons is limited to distances 
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up to 50 lattice parameters (~ 15 /im) in state-of-the-art structures, thereby impeding large-scale 
applications such as integrated circuits. Conversely, scattering in pho 
novel physics such as Anderson localization and non-classical diffusion. 
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The promise of full control over emission and propagation of light has led to a widespread 
pursuit of photonic crystals in recent years [l| . Photonic crystals are dielectric structures in 
which the refractive index varies periodically over length scales comparable to the wavelength 
of light. For three-dimensional periodicities, such crystals promise a photonic band gap, 
i.e., a frequency range for which emission and propagation of light are completely forbidden. 
Ideally, photonic band gap crystals will form a backbone in which many photonic devices, 
such as ultrasmall waveguides, cavities and light sources, are combined to create optical 
integrated circuits [2]. This requires photonic crystals with negligible optical extinction over 
millimeter distances [2]- 

Tremendous progress has been made in the fabrication of photonic bandgap materials 
of the required high refractive index- materials jsl y, |5|, |6| , with low point and plane defect 
densities J6]. Structural variations in size and position of the building blocks, however, are 
intrinsic to three- and two-dimensional (3D, resp. 2D) photonic crystals alike, amounting 
to at least 2 to 7% of the lattice spacing in all current state-of-the-art photonic crystals y, 

n n 

^. While displacements are well-known in condensed matter jj], size polydispersity of 
individual unit cell building blocks, including roughness, is intrinsic to meta-materials such as 
photonic crystals. All such variations can ultimately be traced back to basic thermodynamic 
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but are at present probably limited by materials science. These deviations 
from perfect periodicity cause scattering, and hence exponential attenuation of coherent 
beams propagating through photonic crystals over lengths i, also known as the ' (extinction) 
mean free path'. After propagating over a distance i, a coherent light beam is converted to 
a diffuse glow that corrupts the functionality of any photonic integrated circuit. Conversely, 
short mean free paths open up novel physics related to diffusion of light and ultimately 
Anderson localization of light [9|, |lO|. Therefore, it is crucial to obtain the relation between 
the extinction length i and the structural disorder. In this paper, we derive such a relation 
and test it against available experimental results [11 1. 

We consider extinction in photonic crystals due to scattering by size polydispersity and 
displacements from lattice sites of the structural units (size r) that compose the unit cell 
(lattice spacing a). Light scattering is caused only by the difference in refractive index 
profile of the displaced, slightly polydisperse building blocks as compared to the ideally 
ordered structure. As illustrated in Fig. ^ this difference is a collection of thin shells of 
high and low index material. The polydispersity and displacements of the building blocks 
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FIG. 1: (Schematic) Any 2D or 3D real photonic crystal is an ordered stack of building blocks with 
a spread Ar in their average radius r, each slightly displaced (displacement An) from the lattice 
sites. The real structure is the sum of the perfect crystal and the difference between the real and 
perfect structure. This difference is a collection of thin shells that each scatter weakly. Due to 
their high number density, the shells dominate the scattering loss. 

translate linearly into the shell thickness Ar. Since in many photonic crystals, such as cubic 
(3D) or hexagonal (2D) structures, light transport is isotropic, we treat the ideal crystal as 
an effectively homogeneous medium with index Ues equal to the volume-averaged refractive 
index |l^. Within this framework, the inverse extinction length 
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is the product of three factors [13[: Rayleigh's extinction cross section aRayicigh of each 
shell, the number density of shells p, and a wavelength-dependent geometrical factor F 
which embodies corrections beyond Rayleigh scattering [l^. Since the volume of each shell 
is proportional to its thickness Ar, Rayleigh's extinction cross section is proportional to 
{m — l)^Ar^, where m is the index contrast relative to the background medium. Even 
though scattering by each shell is generally weak, the huge density p set by the number of 
structural units per unit cell causes the scattering mechanism to be important. For Rayleigh 
scatterers, in the low-frequency limit, the dimensionless factor F equ als unity. For weakly 
scattering shells, the Rayleigh-Gans approach is suited to find F [l3|, |l^ . 

We now focus on the extinction length in 3D photonic crystals that consist of spheres 
(mean radius r), such as opals and inverse opals where many data are available. Size 
polydispersity results in scattering due to thin spherical shells with a Gaussian distribution of 
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FIG. 2: (color) Symbols: mean free path [ij] i. in units of a versus normalized frequency a/A 
in polystyrene opals (A) and titania inverse opals (B). Open symbols in (B) were obtained by 
averaging for each a/A total transmission spectra for many samples with different a. The blue 
shaded area indicates the standard deviation. In the stop gap (orange bar), total transmission 
is reduced in excess of l/L due to Bragg reflection of the input beam. This affects the data in 
this limited range. Previous data show that I is unaffected if the frequency is tuned through a 
gap [la]. In both (A) and (B), the extinction length agrees well with the model (|1I2J) with 5R ~ 5% 
(red curves). Green curves represent scaling of £ with w~^, and illustrate the failure of Rayleigh 
scattering models. 

thicknesses. The inverse extinction length £ scales quadratically with the size polydispersity 
Ar and with m — 1, since Rayleigh's extinction cross section for a shell of thickness Ar 
reads crRayieigh = (327r/3)(m — l)^A;ggr^Ar^ (cf. [13j), with k^f^ the wave vector in the effective 
medium. We find that the Rayleigh-Gans correction |l3l | 
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reduces the well-known four th p ower Rayleigh increase of extinction to a nearly quadratic 
dependence on wave vector [15[ . We have checked the validity of our result using the exact 
Mie-solution for spherical shells. Although ioi m> 2 and ke^r > 1 the Rayleigh-Gans result 
underestimates the extinction loss compared to Mie theory, the Mie-model reproduces the 
quadratic scaling with frequency and shell thickness. Our model captures both the effect of 
polydispersity Ar/r and displacements Au/r: calculations of F show that both effects are 
similar in magnitude, and can be combined by taking an effective shell thickness Ar + O.SAw. 
From now on, 6R indicates effective shell thicknesses normalized by the shell radius. An 
essential result of our paper is that given the current fabrication accuracies of 6R ~ 5%, the 



maximuiii extinction length £ is only 50 lattice spacings in high-index crystals at relevant 
frequencies. 

Enhanced backscattering measurements obtained earlier by us have allowed us to deter- 



mine the mean free path d, 1^ in synthetic opals, i.e., fee crystals of close packed polystyrene 
spheres with n = 1.59 and ries = 1-45 [16]. In Figure |21^A), we plot i for a wide normalized 
frequency range, obtained with A = 632, 685 and 780 nm, and many different a. We see 
that i decreases from 100a for frequencies below first order diffraction, to 5a at the highest 
frequencies, where we have converted the wavevector from the scattering model to the fre- 
quency scale a/A typical of photonic crystals. The data and our model agree well on both 
the observed decrease of i with a/A and the magnitude of i, which confirms that extinction 
is due to non-uniformities and displacements of the spheres, assuming 6R = 5%. This value 
matches well with the cumulative effect of polydispersity ~ 2% and rms displacements of 
spheres from their lattice sites (< 3.5% of the nearest neighbor distance), as independently 
determined by small angle X-ray scattering |l^. In contrast, the data refute the often as- 
sumed Rayleigh u'^ dependence ^ ll2|. The degree of extinction is also inconsistent with 
the common assumption that scattering is due to point defects, e.g., missing spheres: From 
the cross-section of a sphere [13| we calculate that the observed scattering would require 
a density of missing spheres larger than 0.13a~^, an order of magnitude larger than the 
estimated density O.Ola"^ p,ll2|- 

We have carried out new experiments to probe scattering losses in photonic crystals with 
high photonic interaction strength, i.e. inverse opals in a Ti02 backbone. The strength of 
the interaction of a photonic crystal with light is gauged by the relative bandwidth S of 
the lowest order gap in the dispersion relation, see Ref. [l|, p. 194. The generally pursued 
large interaction strengths require a large index contrast nhigh/'^-iow and are thus associated 
with stronger scattering, due to the factor {m — 1)^ in Rayleigh's cross section. While the 
magnitude of the non- uniformities is similar to those in the direct opals p| , the inverse opals 
present a much larger index contrast n = 2.7 ± 0.4 (riefr ~ 1.18). We have determined 
the frequency dependence of £ from total diffuse transmission (T = i/L, with L the sample 



thickness 1J|)- We used white-light FTIR spectroscopy to cover a wide normalized frequency 
range for many samples with a = 650 nm to 930 nm. To obtain the absolute magnitudes 
of the mean free paths we calibrated our measurements by measuring the absolute values of 
the transmission (closed symbols) and using enhanced-backscattering data |lq. Figure E^B) 



shows that £ decreases from 100a at a/A = 0.4 to only 4a at a/A = 1.6. This decrease of i. 
for the inverse opals is in excellent correspondence with our prediction (solid curve), taking 
a non-uniformity 5R = 4% that is consistent with independent structural data J3[ . 

To further test the validity of our model, we have analyzed transmission data reported in 
many papers encompassing fee and bcc photonic crystals, with sphere volume fractions from 
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(f = 0.7% to 74% and index contrasts nhigh/''^iow from 1.05 to 
Extinction causes the coherent-beam transmission outside stop gaps to decrease according 
to Lambert-Beer's law Tcoh = e~^^^. In all cases, except of course for the dilute crystal 22 1, 
fitting a power law dependence i oc cj"^ to each data set shows that extinction does not 
increase according to Rayleigh's law. Indeed, Table E] shows that we find exponents x < 4 in 
reasonable correspondence to the exponents predicted by our model in the same frequency 
windows. Similar exponents were recently also observed in Ref. J26|. Fits to our model 
further show that extinction lengths for the wide range of crystals agree with 6R > 4%, 
consistent with typical sphere polydispersities and displacements of 2-5%, as reported in 
Table HI The quantitative agreement of i with Eqs. ()ip2|l confirms that polydispersity and 
displacements of unit cell building blocks determine scattering loss in 3D photonic crystals. 
Given the success of our model, we can now use it to infer the general dependence of 
the extinction length i/a on the photonic interaction strength S and the non-uniformity 
6R. In Figure El we present both i/a and S that are calculated as a function of index 
contrast {m — 1). It is clear that the extinction length decreases both with increasing 



TABLE I: Photonic interaction strength S, structure and extinction in 3D photonic crystals. 



Ref. 
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^sphere / ^intci 


" r/a 


£/a* 


Xc^p{x rg)" 


5R'' 


[22]« 


0.7% 


1.59/1.33 


o.iwf 


105 


4(3.3) 


12% 


[21] "^ 


< 1% 


1.42/1.48 


cp9 


3000 


3(2.6) 


6% 


[2.0] = 


1.6% 


1.59/1.33 


0.143 


1000 


3.3(3) 


15% 


[17] <= 


2% 


1.32/1.47 


cp 


700 


2.6(2.6) 


6% 


[24]' 


3% 


1.59/1.33 


cp 


100 


< 2(2.5) 


7% 


[251" 


5% 


1.41/1.0 


cp 


17 




9% 


[23]" 


5.5% 


1.45/1.0 


cp 




2.6(2.4) 




Fig.HIa)'' 


7% 


1.59/1.0 


cp 


50 


1.8(2.4) 


5% 


Fig.lSb)' 


11% 


1.0/2.7 


cp 


40 


2.6(2.5) 


4% 



"'^sphere /™intor: refractive indices of spheres, resp. background medium. 'The £/a are for a/A in the first stop gap. '^Decay 

powers X obtained by fitting £ oc lij~^ to the data (xcxp), resp. model Eq. I1I2I in the same frequency range {xhq). ''Effective 

shell radii that best fit the data over the full available frequency range. "^Transmission. ■'Bcc instead of fee. 

9cp=close-packing, r/a = l/\/8. ''Enhanced backscattering. 'Diffuse total transmission. 



photonic strength and with increasing structural disorder. We also present the experimental 
extinction data from Table Ufor fee opals and inverse opals, showing again a good agreement 
with our model with 5R > 4%. A photonic band gap requires interaction strengths beyond 
S = 0.15; extinction lengths less than 20 lattice spacings are expected at the current level of 
fabrication accuracy. Ultimately, one hopes to realize photonic crystals that combine many 
optical functions. Recent technology roadmaps foresee crystals containing ~ 10^ optical 
functions per mm^ (Ref. |2|, P- 245), requiring negligible loss over more than mm distances. 
From the general scaling of extinction with non-uniformity we conclude that applications 
of photonic band gap crystals in circuits require a formidable tenfold increased perfection 
in statistical fabrication accuracy to 6R < 0.25%, or subnanometer precision. Such an 
improvement is far beyond the current state-of-the-art [jj, |2| • 
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FIG. 3: (color) Universal dependence of extinction length on interaction strength S in photonic 
crystals. Dashed curves: extinction calculated for fee inverse opals (assuming 26% high index 
material). Symbols: our data (a), and literature data analyzed by us (♦). Observed losses are 
consistent with 5R > 4%. If i is shorter than the length needed for Bragg diffraction, structures are 
essentially disordered (shaded red). Complete band gaps are expected for S > 15% (shaded green). 
Photonic crystal integrated circuits require i/a > 10^ at 5 > 15%, beyond current state-of-the-art. 



Although 3D photonic crystals potentially offer the best platform for photonic crystal 
functionality, 2D photonic crystals possess many of the desired properties with the advantage 
of ease of fabrication. While the fabrication methods are radically different, 2D photonic 
crystals suffer from similar polydispersity and displacements of their unit cell building blocks 
as 3D crystals pi] . To obtain the scattering losses, we consider 2D crystals of infinitely long 
cylinders. Now, Rayleigh's cross section per unit length crRayicigh = (37r^)(m — lyk^^r'^Ar'^ 
of thin cylindrical shells of thickness Ar and radius r increases with the cube of the optical 
frequency [l3| . In the relevant range of cylinder radii, the Rayleigh-Gans model causes the 
ui~^ dependence of i in the Rayleigh-limit to be reduced to a;~^'^ since |l^ 

F{hsr) ^ 0.488(A;effr)-o-^ (3) 



For a hexagona 
silicon crystals 



lattice of air cylinders in silicon with r/a = 0.45, typical for macroporous 
27[, we find i ~ 40a for frequencies near lowest order stop gaps, assuming a 
non-uniformity 6R of 5%. A much larger i is required for integrated circuit applications. 

Many efforts currently focus on quantifying losses in 2D crystals made from high index 
slabs on lower index cladding layers, for which the nonuniformity SR is around 5% [l|, |8| . 
Although the guided wave profile normal to the slab is not incorporated in our model, we 
believe that Eq. yields a reasonable estimate of scattering due to nonuniformity of the 
air holes in such structures. Similar to 3D, applications of 2D structures in photonic crystal 
integrated circuits require a formidable increase in fabrication accuracies beyond the current 
state-of-the-art [ij, l2] • These scattering losses add to currently widely studied out-of-p ane 
scattering that is intrinsic even to hypothetical perfectly fabricated 2D crystal designs 28|. 
In contrast to out-of-plane loss, however, statistical variations cannot be reduced by design 
optimization. 

Scattering in photonic crystals opens opportunities to explore new phenomena in multiple 
scattering of light [10|. Photonic crystals allow unique control over fundamental aspects, 
such as the transport velocity or anisotropics of light diffusion. A fascinating application 
is the possibility to localize light, which could serve to enhance non-linear interactions J9j. 
According to the usual loffe- Kegel criterion, Anderson localization occurs when the mean free 
path is so strongly reduced that its product with the wave vector equals one: A;£ ~ 1. It has 
been proposed that in photonic crystals this challenging criterion is relaxed to ki ~ 1/^ 
with ps the modification of the photonic density of states (DOS) relative to free space J2£ 

8 



Since the DOS is strongly reduced in photonic gaps, locahzation of hght may even be feasible 
with the relatively long mean free paths predicted by our model. 
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